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Abstract. It is shown how to construct exactly gauge-invariant S'-matrix elements 
for processes involving unstable gauge particles such as the Z° boson. The results are 
applied to derive a physically meaningful expression for the cross-section a(e + e~ —> 
Z°Z° ) and thereby provide a solution to the long-standing problem of the unstable 
particle. 


I INTRODUCTION 

A resonance is fundamentally a non-perturbative object and is thus not amenable 
to the methods of standard perturbation theory. In order to describe physics at 
the Z° resonance one is forced, therefore, to employ some sort of non-perturbative 
procedure. Such a procedure is Dyson summation that sums strings of one-particle 
irreducible (1PI) self-energy diagrams as a geometric series to all orders in the 
coupling constant, a and effectively replaces the tree-level Z° propagator by a 
dressed propagator, 

i ^ i v ( n zz {s) \ n = i 

s — Mz s — M'z n \s — Mz) s — Mf — II^(s) 

where Uz Z (s) is the one-loop Z° self-energy. The problem here is that electroweak 
physics is described by a gauge theory. Results of calculations of physical processes 
must be exactly gauge-invariant but this comes about through delicate cancella¬ 
tions between many different Feynman diagrams each of which is separately gauge- 
dependent. The cancellation happens at each order in a when all diagrams of a 
given order are combined. The Z° self-energy, I Jzz(s), is gauge-dependent at O(a) 
and hence the rhs of eq.(l) is gauge-dependent at all orders in a. 

If the dressed propagator is used in a finite-order calculation the result will be 
gauge-dependent at some order because the will be no diagrams available to cancel 






the gauge-dependence beyond the order being calculated. This gauge-dependence 
should be viewed as an indicator that the approximation scheme being used is 
inconsistent or does not represent a physical observable. 

In constructing amplitudes that represent physical observables, one must take 
care to respect the requirements that are laid down by analytic A-matrix theory 
[1], These conditions are derived from general considerations such as energy con¬ 
servation and causality and it will be found that appealing to them leads to a 
procedure for generating gauge-invariant amplitudes with no flexibility in the final 
result. The procedure will also make it possible to give an answer to the problem of 
the unstable particle that was put forward by Peierls in the early fifties [2], Much 
of what appears here can be found in ref.s [3-5]. 


II GAUGE-INVARIANT 5-MATRIX ELEMENTS 

We will begin by reviewing what is known about S'-matrix elements for processes 
containing unstable particles. An unstable particle is associated with a pole, s p , 
lying on the second Riemann sheet below the real s axis. The scattering amplitude, 
A(s) for a process which contains an intermediate unstable particle can then be 
written in the form, 


A(s) =-f-+ B(s) (2) 

where R and s p are complex constants and B(s) is regular at s = s p . The first 
term on the rhs of eq.(2) will be called the resonant term and the second is the 
non-resonant background term. It is known that s p is process-independent in the 
sense that any process that contains a given unstable particle as an intermediate 
state will have its pole at the same position. It can also be shown from Fredholm 
theory that the residue factorizes as, R — Rj-Rf, into pieces that depend separately 
on the initial- and final-state. One can prove, by very simple arguments, that s p , 
R and B(s ) are separately and exactly gauge-invariant [3]. 

^Frorn analytic S'-matrix theory it is known that production thresholds are asso¬ 
ciated with a branch cut. Branch points for stable particles lie on the real s-axis and 
those for unstable particles, such as W bosons, lie significantly below it. Provided 
there are no nearby thresholds the amplitude, A(s), can be adequately described 
in the resonance region by a Laurent expansion about the pole, s p . It should be 
emphasized that A(s) can always be written in the form (2) even when thresholds 
are present. In that case B(s) will be an analytic function containing a branch 
point and the rhs of (2) continues to be an exact representation of A(s). Laurent 
expansion refers to how the resonant term is identified. The aim is only to separate 
A(s) into gauge-invariant resonant and background pieces and there is no necessity 
to perform a Laurent expansion beyond its leading term although this can provide 
a useful way for parameterizing electroweak data [6]. 



Let us consider the production process for a massless fermion pair in e + e _ —> //. 
Away from the Z° resonance it is known how to calculate the amplitude to arbitrary 
accuracy in a gauge-invariant manner. Near resonance we must perform a Dyson 
summation and then separate the amplitude into its resonant and background 
pieces. Doing so allows expressions for the gauge-invariant quantities, s p , R and 
B(s) to be identified in terms of the 1PI functions that occur in perturbation theory. 
The scattering amplitude for e + e _ —» // is 
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The pole position s p is a solution of the equation 
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and F Z z(s) is defined through the relation 
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It should be emphasized that eq.(3) is exact and valid anywhere on the complex 
s-plane. The effect of Z °-7 mixing has been included. The quantity II 77 (s) and 
n Z 7 (s) are the photon self-energy and the Z- 7 mixing respectively. Viz(q 2 ), Vzf{q 2 ) 
are the initial- and final-state Z° vertices, into which the external wavefunctions 
have been absorbed, and V^(q 2 ), V^f(q 2 ) are the corresponding vertices for the 
photon. Here B(s,t) denotes 1PI corrections to the matrix element that include 
things like as box diagrams. The first term on the rhs of eq.(3) is the resonant 
part of A(s, t) and the three terms on the second line taken together are form the 
non-resonant background. 

Calculations are most conveniently performed in terms of the real renormalized 
parameters of the theory such as the renormalized mass, M Z - Eq.( 6 ) can be solved 
iteratively in terms of Mz to give 


s p — M% + Bzz(M'z) + ... 


( 8 ) 



The rhs of eq.(8) may be substituted for s p where it appears in eq.(3)-(5). Taylor 
series expansion can then be used to obtain perturbative expressions for s p , R and 
B (s) in terms of Greens functions with real arguments up to any desired order. At 
any given order these perturbative expressions will be exactly gauge-invariant as 
will scattering amplitudes constructed from them. 

A couple of points should be noted here. The appearance of Greens functions 
with complex arguments in eq.(3) is a natural consequence of the analyticity of the 
S-matrix. The A-matrix itself is never evaluated with complex s. The arguments 
of A(s, t) on the lhs of eq.(3) is real as is the ‘s’ in the denominator of the first term 
on the rhs. 

In the procedure described above, one starts by extracting the resonant term 
in a scattering amplitude by Laurent expansion about the exact pole position s p 
and then specializes to lower orders by further expanding about the renormalized 
mass. Other authors [7,8] have attempted to apply the techniques described above 
by first expanded about the renormalized mass and then added a finite width in 
the denominator of the resonant part by hand. That procedure cannot be justified 
and leads to problems when one treats processes like e + e“ —> W + W~. It gives rise 
to spurious threshold singularities or complex scattering angles due the production 
threshold’s being incorrectly located on the real axis. In section IV the process 
e + e _ —»• Z°Z° will be treated and no threshold singularities or complex scattering 
angles will arise. 


Ill THE PROBLEM OF THE UNSTABLE PARTICLE 


We have thus succeeded in our goal of producing exactly gauge-invariant scat¬ 
tering amplitudes to arbitrary order. One might ask whether what has been done 
is just a mathematical trick, in which case the gauge-invariance is accidental, or 
does it have some physical interpretation. In this section it will be shown that the 
latter is true. 

Recall that the coordinate space dressed propagator for a scalar particle has an 
integral representation 
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The integrand has a pole at k 2 = s p where s p is a solution of the equation s — 
m 2 + II(s) = 0 and as in the previous section we define F(s) via the relation 
s — m 2 T II(s) = (s — s p )/F(s). The dressed propagator can then be written as 
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that separates resonant and non-resonant pieces. Performing the ko integration 
resonant gives 
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where A; 0 = yk 2 + s p . The non-resonant term contributes only for t — t' and 
so represents a contact interaction. The resonant part spits into two terms that 
contribute when t > t' or t < t' and therefore connects points x and x' that are 
separated by a finite distance in space-time. 

The problem of the unstable particle [2] is may be roughly stated as follows: 5- 
matrix theory deals with asymptotic in-states and out-states that propagate from 
and to infinity. Unstable particles cannot exist as asymptotic states because they 
decay a finite distance from the interaction region. Indeed it is known [9] that the 
5-matrix is unitary on the Hilbert space spanned by stable particle states and hence 
there is not even any room to accommodate unstable particles as external states. 
How can one use the 5-matrix to calculate, say, the production cross-section for an 
unstable particle when it cannot exist as an asymptotic state? 

In the first part of this section it was shown that the resonant part of the dressed 
propagator connected points with a finite space-time separation. When a similar 
analysis is applied to a physical matrix element, such as eq.(3), one concludes that 
the resonant part describes a process in which there is a finite space-time separation 
between the initial-state vertex, U, and the final-state vertex, Vf. In other words, 
the resonant term describes the finite propagation of a physical Z° boson. The non¬ 
resonant background represents prompt production of the final state. As these two 
possibilities are, in principle, physically distinguishable, they must be separately 
gauge-invariant. 

We can thus use finite propagation as a tag for identifying unstable particles 
without requiring that they appear in the final state. This is, after all, the way b- 
quarks are identified in vertex detectors. A production cross-section for an unstable 
particle is obtained by extracting the resonant part of the matrix element for a 
process containing that particle in an intermediate state and summing over all 
possible decay modes. 


IV THE PROCESS e+e -> Z°Z° 

In this section we will calculate the cross-section for e + e~ —> Z°Z°. This is of 
both theoretical and practical importance. On the theoretical side it represents 
an example of a calculation of the production cross-section for unstable particles. 
On the practical side, at high energies e + e _ —»• Z°Z° will be a dominant source of 
fermion pairs (/i/i) and (f 2 f 2 ) due to its double resonant enhancement and hence 



i x(e + e~ —> Z°Z °) is an excellent approximation to the cross-section for 4-fermion 
pair production. If the experimental situation warrants it, background terms can 
also be included without difficulty. 

In the case of e + e _ —> //, dealt with in section II, the invariant mass squared of 
the Z°, s, is fixed by the momenta of the incoming e + e _ . For the process e + e~ —» 
Z°Z° the invariant mass of the produced Z 0, s is not constant and must be somehow 
included in phase space integrations. It is not immediately clear how to do this and 
without further guidance from S'-matrix theory there would seem to be considerable 
flexibility in how to proceed. A new ingredient is required and that is to realize that 
an expression for an S-matrix element can always be divided into a part that is a 
Lorentz-invariant function of the kinematic invariants of the problem and Lorentz- 
covariant objects, such as j) etc. The latter are known as standard covariants 
[10-13]. It is the Lorentz invariant part that satisfies the requirements of analytic 
S'-matrix theory and from which the resonant and non-resonant background parts 
are extracted while the Lorentz covariant part is untouched. 

To calculate a(e + e~ —> Z°Z°), we begin by constructing the cross-section for 
e + e _ —> Z°Z° —> (/i/i)(/ 2 / 2 ), and will eventually sum over all fermion species. The 
part of the full matrix element that can give rise to doubly resonant contributions 
can be written as 

M = Y^lVe+T^Ue-jMifa^fap l) 
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where T^ JU are Lorentz covariant tensors that span the tensor structure of the ma¬ 
trix element and y^, 7 r are the usual helicity projection operators. The squared 
invariant masses of the f\ f\ and / 2 /2 pairs are p\ and p\. The M t , II^ and Vzf are 
Lorentz scalars that are analytic functions of the independent kinematic Lorentz 
invariants of the problem. 

To extract the piece of the matrix element that corresponds to finite propagation 
of both Z 0, s we extract the leading term in a Laurent expansion in p\ and p\ of 
the analytic Lorentz-invariant part of eq.(12) leaving the Lorentz-covariant part 
untouched. This is the doubly-resonant term and is given by 

.AA y ^ \v e +T lw u e -]Mi(t, u, s p , s.p) 
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where F Z z defined by a relation like (7). It should be emphasized that eq.(13) is 
the exact form of the doubly-resonant matrix element to all orders in perturbation 



theory that we will now specialize to leading order. It is free of threshold singular¬ 
ities noted that were found by other authors [ 8 ]. In lowest order eq.(13) becomes, 
up to overall multiplicative factors, 
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where = 7 ^(^ e - - /Ji) 7 „, Af x = t \ T 2 V = 7 v (j) e - - j4 2 )7 m , M 2 = u 1 and the 
final state vertex corrections take the form Vzf L = iefd^pL and Vzf R = iePn'yR, The 
left- and right-handed couplings of the Z° to a fermion / are 
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Squaring the matrix element and integrating over the final state momenta for 
fixed p\ and p\ gives 
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where the sum is over fermion species. Note that p(p 2 ) —■► 5(p 2 — M^)0(po) as 
Irn(,Sp) —» 0 which is the result obtained by cutting a free propagator. The variables 
s, t, u, pi and p 2 in eq.(15) arise from products of standard covariants and external 
wave functions and therefore take real values dictated by the kinematics. 

Integrating over t, pi and pi leads to 
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where 
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and A = + pf + ~~ 2sp 2 — 2sp| — ^PiPl- For p\ = p\ = M| this agrees with 

known results [14]. 
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